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ABSTRACT 


Similarity  oe  signal  waveform  across  the  Large  Aperture 
Seismic  Array  (LASA)  was  studied.  The  analysis  technique  depended  on 
differences  in  waveform  shape  but  not  on  amplitude  differences. 

The  waveform  was  found  to  be  very  similar  both  within 
subarrays  and,  except  for  a  few  cases,  between  subarrays.  Thus,  1-pt 
(amplitude)  equalization  usually  is  sufficient  when  proces  ng  LASA  data 
both  on  the  subarray  and  large-array  levels. 
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SECTION  I 
INTRODUCTION 


The  purpose  of  this  study  was  to  measure  the  similarity  of 
signal  waveforms  across  LASA  both  within  subarrays  and  between  subarrays. 
The  following  three  types  of  signal  similarity  were  analyzed: 

(1)  Single  seismometers  within  a  subarray.  All  seismo¬ 
meters  at  each  of  two  subarrays  (F3  and  F4)  were 
used  for  this  study. 

(2)  Single  seismometers  between  subarrays.  The  center 
seismometers  from  all  subarrays  were  used  for  this 
study. 

(3)  Subarray  outputs.  The  time -shifted  sums  from  all 
subarraye  were  used  for  this  study. 

The  five  events  used  for  analyses  (1)  and  (2)  are  listed  in 
Table  1;  the  20  events  used  for  analysis  (3)  are  listed  in  Table  2.  Two 
events  are  common  to  both  tables.  For  analyses  (1)  and  (2),  the  events  were 
chosen  from  a  suite  which  had  been  demultiplexed  under  another  task  of  the 
LASA  contract.  For  analysis  (3),  the  events  were  chosen  from  the  suite 
used  to  study  the  relative  capabilities  of  large  and  small  seismic  arrays 
for  event  identification. 

All  events  were  resampled  to  a  0. 1  -sec  rate  and  bandpass 
filtered  with  a  zero-phase  0.  8  to  2.  8  cps  digital  filter  (Figure  1)  to  reduce 
the  low-frequency  ambient  noise.  Only  events  with  a  large  signal -to -noise 
ratio  were  used. 


Texas  Instruments  Incorporated,  1967:  Large  Array  Signal  and  Noise 
Analysis,  A  study  of  the  Relative  Capabilitv  of  Large  and  Small  Seismic 
Arrays  for  Event  Identification,  Spec.  Rpt.  No.  1,  Contract  AF33(657) -16678, 
20  April. 
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Amplitude  Response  of  3.7-sec 
Zero-Phase  Bandpass  Filter: 
100  msec  Sample  Rate;  0.  8  to 
2.  8  cps  Passband 


SECTION  II 

ANALYSIS  TECHNIQUE 


Waveform  similarity  was  measured  by  computing  correlation 
coefficients  between  given  traces  anl  a  reference  trace.  For  single  seis¬ 
mometers  within  a  subarray,  the  reference  trace  used  was  the  subarrav 
time-shift-and-sum.  The  time-shift-and-sum  of  the  single  seismometers 
involved  was  the  reference  trace  used  for  single  seismometers  between 
subarrays.  For  subarray  outputs,  the  reference  trace  used  was  the  LASA 
time-shift-and-sum.  In  each  case,  the  reference  trace  was  thus  the  appro¬ 
priate  average  of  the  input  traces. 

Using  the  appropriate  average  as  a  reference  rather  than  some 
arbitrary  individual  trace  was  preferable  because,  in  cases  where  that  indi¬ 
vidual  trace  had  a  significantly  different  waveform  than  the  other  traces,  a 
low  set  of  correlation  coefficients  resulted.  By  using  the  average  trace  as 
?  reference,  only  ne  correlation  coefficient  associated  with  that  trace  was 
low.  The  correlation  coefficient  p  was  defined  as 

cp.r(T) 

p  *  —  - -  - 

V^a'01  Vv,o) 

where 

<Pir(T)  is  the  maximum  lag  of  the  cu  ss correlation  between 
the  individual  and  reference  traces 

tPu(0).  ^rr(°)  is  the  zero-lag  autocorrelations  of  the  individual 
and  reference  traces 


The  "full-house"  correlation  technique  was  used  in  the  computations.  That 
is,  the  gateti  were  chosen  in  such  a  way  that  p  =  1  if  the  waveforms  are 
identical  except  for  being  displaced  by  a  lag  T. 

The  gate  lengths  used  to  compute  the  correlations  were  visually 
chosen  to  include  the  main  signal  arrival.  When  the  signal  had  fairly  long 
duration,  several  gates  were  usually  selected. 

The  correlation  coefficient  was  chosen  as  the  analysis  tool 
because  it  measured  duierences  in  waveform,  was  a  relatively  inexpensive 
method  of  analysis,  and  did  not  depend  on  g  lin  differences  between  channels. 
However,  the  effect  of  a  linear  filter  was  reflected  in  the  correlation  co¬ 
efficient  (as  contrasted  with  the  2-channel  coherence  technique).  That  is, 
let 


i  (t)  =  r  (s)  *  h  (s) 


where 


Then, 


Mt) 

is  an  individual  trace 

r(s) 

is  the  reference  trace 

h(s) 

is  a  linear  fil  j;* 

* 

stands  for  convolution 

CD.  It) 
ir 

a  r  (s)  §  [r  (s)  *  h  (s)| 

a  Cp  (t)  *  h  (t) 
rr 

<Pu(t) 

a  fr{s)  *  h{s)j  #  [r(s)  * 

*  v^’V" 


CP..{0)  *  Cp  (T)  *  Cp.  ,{T) 

u  rr  hh 


T  a  0 


where  ♦  stands  for  "correlated  with.  " 
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cp  (t)  *  h(t) 

_ _ _ rr _ 

VV<°>  ^ rr(l)*  V(t> 


and,  unless  h(s)  is  a  1-pt  filter,  p  ^  1. 

The  mean  and  variance  of  the  set  of  correlation  coefficients 
were  computed  for  each  event,  and  measures  of  both  signal  degradation  and 
signal -to -noise  ratio  were  obtained.  Signal  degradation  "L"  was  obtained  by- 
choosing  the  largest  peak-tto-peak  amplitude  in  the  first  few  cycles  and 
computing 


where 


L  *  20  log1() 


A. 

l 


i  =  1 


L  is  the  signal  degradation  in  db 

is  the  peak-to-peak  amplitude 
on  the  reference  trace 

is  the  peak-to-peak  amplitude  on  an 
individual  trace 

N  is  the  number  of  individual  traces 


_s_ 

N 


max 


N 


rms 


where 


Amax  iS  the  maximum  zero-to-peak  amplitude 
in  the  first  few  cycles 


N 


nus  is  tlie  rms  level  of  the  noise  immediately 
preceding  the  signal 


Finally,  an  attempt  was  made  to  relate  the  means  and/or 
variances  of  the  correlation  coefficients  to  both  signal  degradation  and 
signal-to-noise  ratio. 
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SECTION  III 

PRESENTATION  OF  DATA 


A.  INTRASUBARRAY 


To  analyze  the  similarity  of  single -seismometer  outputs 
within  a  subarray,  the  five  events  listed  in  Table  1  were  used. 

Subarrays  F3  and  F4  were  chosen  because  they  usually  had  larger 
signal -to -noise  ratios  than  other  subarrays.  Figure  2  shows  the 
reference  trace  {subarray  time -shift -and-sum)  for  each  event  and 
the  gates  used  in  computing  the  correlation  coefficients. 

Table  3  gives  the  correlation  coefficients  of  each 
seismometer  for  each  event.  It  can  be  seen  that  most  of  the  coefficients 
were  close  to  1.0.  Excluding  subarray  F3  for  the  Andreanof  Islands 
event,  only  seven  coefficients  were  less  than  0.  9  and  only  one  (seismom¬ 
eter  45,  subarray  F3,  Colombia  event)  was  less  than  0.  80.  Subarray  F3 
was  anomalous  for  the  Andreanof  Islands  event. 

Figure  2  shows  that  the  large -amplitude  arrival  for 
subarray  F3  was  delayed.  It  was  the  first  arrival  for  the  other 
subarrays.  Because  of  this  delay,  the  correlation  coefficient  for  F3 
was  computed  over  a  gate  with  a  relatively  small  signal -to-noise 
ratio,  adversely  affecting  the  correlation  coefficient.  A  second  set  of 
coefficients  computed  over  a  gate  which  contained  the  large  arrival 
was  considerably  higher.  Figure  3  shows  thf-  Panama  event  as 
recorded  at  both  subarrays  and  is  typical  of  the  waveform  duplication 
observed  at  the  subarray  level. 

The  maximum  variations  in  amplitude  (i.  e. ,  the  ratio  of 
the  largest  to  smallest  amplitude)  across  a  subarray  for  each  event  are 
given  in  Table  4.  Variations  ranged  up  to  7  db,  showing  a  need  for 
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Tablj  3 

CORRELATION  COEFFICIENTS  WITHIN  SUBARRAYS 


Seismometer 

Andreanof 

Islands 

Colombia 

I 

Mexico 

Pakistan 

Panama 

10 

0.976 

0.977 

0.976 

0.977 

0.974 

21 

0.979 

0.  904 

0.993 

0.945 

0.997 

41 

0.  885 

0.978 

0.986 

0.994 

0.963 

61 

0.855 

0.976 

0.980 

0.972 

0.964 

81 

0.  868 

0.866 

0.958 

0.958 

0.980 

22 

0.945 

0.978 

0.982 

0.998 

0.944 

32 

0.853 

0.972 

0.993 

0.943 

0.987 

52 

0.922 

0.976 

0.948 

0.957 

0.994 

72 

0.937 

0.949 

0.982 

0.982 

0.980 

23 

0.914 

0.977 

0.978 

0.994 

0.984 

43 

0.947 

0.955 

0.964 

0.998 

0.951 

63 

O' 

o 

00 

• 

o 

0.893 

0.994 

0.992 

0.988 

83 

0.868 

0.903 

0.963 

0.987 

0.997 

24 

0.969 

0.990 

0.995 

0.978 

0.995 

34 

0.921 

0.931 

0.988 

0.997 

9.992 

54 

0.817 

0.971 

0,993 

0.946 

0.983 

74 

0.747 

0.885 

0.984 

0.949 

0.980 

25 

0.978 

0.956 

0.  964 

0.954 

0.968 

45 

0.856 

0.696 

0.985 

0.976 

0.953 

65 

0.887 

0.C71 

0.  V73 

0.959 

0.975 

85 

0.  878 

0.92V 

0.984 

0.989 

0.964 

26 

0.972 

0.94C 

0.991 

0.962 

0.992 

36 

0.865 

0.978 

0.986 

0.952 

0.991 

56 

0.926 

0.893 

0.994 

0.989 

0.982 

76 

0.804 

0.913 

0.940 

0.951 

0.947 

™  LtJS: 
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Table  3  (Contd) 


F4 


Seismometer 

. . . 1 

Andreanof 

Islands 

Colombia 

Mexico 

Pakistan 

Panama 

10 

0.958 

0.945 

0.980 

0.994 

0.994 

21 

0.940 

0.911 

0.957 

0.980 

0.955 

31 

0.956 

0.974 

0  972 

0.992 

0.954 

51 

0.986 

0.967 

0.989 

0.958 

0.990 

71 

0.970 

0.947 

0.961 

0.937 

0.  980 

22 

0.929 

0.937 

0.968 

0.991 

0.995 

42 

0.948 

0.988 

. 

0.980 

0.976 

62 

0.962 

0.954 

0.989 

0.997 

82 

0.992 

0.971 

0.981 

0.995 

0.989 

23 

0.933 

0.994 

0.992 

0.994 

0.973 

33 

0.994 

0.970 

0.  988 

0.997 

0.956 

53 

0.983 

0.972 

0.955 

0.995 

0.996 

73 

0.965 

0.947 

0.951 

0.982 

0.977 

24 

0.954 

0.988 

0.999 

0.993 

0.992 

44 

0.988 

0.948 

0.996 

0.987 

0.989 

64 

0.  974 

0.900 

0.987 

0.957 

0.991 

84 

0.974 

0.925 

0.  991 

0.969 

0.991 

25 

0.988 

0.967 

0.990 

0.997 

0.984 

35 

0.992 

0.933 

0.998 

0.993 

0.998 

55 

0.988 

0.  932 

0.987 

0.963 

0.993 

75 

0.976 

0.908 

0.973 

0„  990 

0.988 

26 

0.994 

0.929 

0.974 

0.991 

0.988 

46 

0.956 

0.942 

0.994 

0.995 

0.994 

66 

0.945 

0.966 

0.961 

0.985 

0.969 

86 

0.988 

0.874 

0.970 

0.973 

0.989 

- - 
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Table  4 

MAXIMUM  AMPLITUDE  VARIATION  ACROSS  SUBARRAYS 


1 - 

Event 

Maximum  Amp 
Variation  (db) 

Maximum 

Seismometer 

Minimum 

Seismometer 

Andreanof  Islands 

7.0 

61 

54 

Colombia 

4.2 

76 

72 

Mexico 

1.7 

81 

72 

c* 

Ui 

W.  Pakistan 

2.6 

34 

76 

Panama 

4.3 

26 

72 

Average 

4-° 

Andreanof  Islands 

6.4 

24 

86 

Colombia 

6.9 

73 

71 

Mexico 

2.9 

84 

86 

* 

Ui 

W.  Pakistan 

3.4 

24 

26 

Panama 

4.7 

24 

86 

Average 

4.9 

0 


amplitude  equalization  at  the  subarray  level.  Note  that  seismometers  with 
the  largest  and  smallest  amplitudes  were  different  for  different  events  — 
possibly  due  to  either  statistical  seismometer -gain  fluctuations  or  a  "tuning" 
of  individual  seismometers  to  epicentral  regions  which  was  probably  caused 
by  nonhomogenous  seismometer -ground  coupling. 

Tabla  5  lists  the  correlation  coefficients  for  the  Pakistan  event 
(Figure  4)  for  several  gate  lengths  and  shows  that  they  remained  high  even  for 
long  gates.  This  indicates  that  scattered  energy  is  not  a  problem  for  this 
event.  Similar  results  for  other  events  indicate  that,  in  general,  scattered 
energy  is  not  a  problem  at  the  LASA  site. 

Table  6  lists  the  means  and  variances  of  the  correlation  coef  - 
ficienti,  the  signal  degradations,  and  the  signal -to -noise  ratios  of  the 
subarray  sums  for  the  five  events. 

Subarray-F4  signals  seem  to  have  more  similarity  than  those 
in  F3,  as  indicated  by  their  higher  means,  lower  variances,  and  smaller 
degradation  values.  Those  events  with  the  larger  coefficient  means  had 
less  signal  degradation  (as  would  be  expected),  although  all  events  had  less 
than  1-db  degradation  .  With  the  exception  of.  the  Andreanof  event  on  F3,  all 
events  had  sufficiently  large  signal -to -noise  ratios  to  prevent  the  ambient 
noise  from  affecting  the  correlation-coefficient  values. 

B.  SINGLE  SEISMOMETERS  BETWEEN  STJ  BARR  AYS 

To  analyze  the  similarity  of  single  seismometer  outputs 
between  subarrays,  the  same  five  events  listed  in  Table  I  were  used. 
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The  center  seismometers  from  all  subarrays  were  chosen,  with 
the  time -shift -and -sum  used  as  the  reference  trace.  Figure  5  shows  the 

reference  trace  for  each  event  and  the  gates  used  in  computing  the  correlation 
coefficients. 

Table  7  gives  the  correlation  coefficient  for  each  seismometer 
for  each  event.  Note  that  the  coefficients  were  generally  large  (greater  than 
0.8),  although  they  were  smaller  on  the  average  than  the  intrasubarray  values. 
However,  a  few  seismometers  had  very  low  values  (e.g. ,  the  D2  and  the  El 

seismometers  for  the  Colombia  event);  the  significance  of  these  low  values 
is  discussed  in  subsection  C. 

Figure  6  shows  the  Mexico  event  as  recorded  by  the  21  center 
seismometers  and  indicates  that  the  waveform  duplication  of  single  seismom¬ 
eters  .across  LASA  was  quite  good.  Maximum  variation  in  amplitude  across 
LASA  naturally  was  much  larger  than  across  a  subarray  (up  to  18  db  as 
compared  to  7  db  for  the  subarrays).  Again,  amplitude  equalization  was 
necessary  (  and,  in  most  cases,  probably  sufficient). 

Table  b  lists  the  correlation -coefficient  means  and  variances, 
the  signal  degradation,  and  the  signaX-to-noise  ratios  on  the  reference  traces 
for  the  five  events.  Means  were  about  0.  08  lower  than  intra-array  means, 
and  variances  were  slightly  higher.  Again,  events  with  larger  coefficient 
means  had  less  signal  degradation. 

C.  SUBARRAY  OUTPUTS 

The  twenty  vents  listed  in  Table  2  were  used  to  analyze  the 
similarity  of  eubarray  ov  -uts.  Figure  7  shows  the  LASA  sum  for  each 
event  and  the  gates  used  in  computing  the  correlation  coefficients. 

Table  9  lists  the  correlation  coefficients  for  each  subarray 
output  for  each  event.  Figure  8  shows  the  location  of  the  events  on  a  polar 
plot  centered  at  LASA.  The  coefficients  were  generally  large  (greater  than 

0.8),  although  smaller  on  the  average  than  the  intrasubarray  values. 
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Table  7 

CORRELATION  COEFFICIENT  FOR 
SINGLE  SEISMOMETERS  BETWEEN  SUBARRAYS 


Aleutian 

0.975 

0.929 

0.776 


0.878 


0.919 


0.901 


0.860 


0.979 


0.945 


0.944 


0.938 


0.933 


0.861 


0.975 


0.983 


0.850 


0.970 


0.933 


0.926 

0,929 


Colombia 

0.855 


0.783 


0.786 


0.798 


0.  877 
0.811 
0.964 
0.936 
0.939 
0.  828 
0.  875 
0.862 


0.760 


0.  844 


0.  338 


0.  875 


0.  760 


0.421 


0.756 
0.  878 


Mexico 


0.771 


0.980 


0.908 


0.923 


0.899 


0.901 


0.910 


0.925 


0.643 


0.  933 


0.970 


0.935 


0.927 


0.916 


0.932 

0.919 


0.931 


0.953 


0.925 


0.  551 


Pakistan 


0.  934 


0.876 


0.971 

0.948 

0.964 


0.857 


0.957 

0.736 

0.941 

0.953 

0.  885 

0.766 

0.929 

0.948 


0.902 


0.877 


0.801 


0.811 

0.877 


Panama 


0.858 


0.938 


0.967 


0.766 


0.943 


0.964 


0.982 


0.994 


0.941 


0.908 


0.773 


0.920 


0.926 


9.973 


0.960 


0.855 


0.974 


0.971 


0.944 

0.  98i 


0.836 


0.941 


0.941 


0.847 


0.926 
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Event  27  (Tadzhik),  which  had  excellent  subarray  waveform 
duplication  across  LASA,  is  shown  in  Figure  9. 

Table  9  shows  some  subarrays  having  very  low  coefficients 
for  a  few  evehts.  For  example,  the  coefficients  for  subarrays  B3  and  D1 
(and,  to  a  lesser  extent,  AO  and  F4)  were  low  for  all  five  nuclear  explosions 
from  Kazahk. 

Figure  10  shows  event  105;  the  signal  waveform  was  significantly 
different  on  these  subarrays. 


Figure  11  compares  the  Di  output  and  LASA  sum  for  all  five 
events.  While  Dl  was  consistently  different  from  the  sum,  it  was  similar  for 
each  event. 


Among  other  anomalous  outputs  were  E3  (very  low  for  event  2, 
Fiji),  and  D3  and  E2  (very  low  for  Event  101,  Algerian  nuclear  blast).  Note 
that  these  subarrays  were  "normal"  for  the  ether  events. 

The  anomalies  seemed  to  depend  quite  critically  on  both  event 
azimuth  and  epiceatral  distance  (i.  e. ,  angle  of  incidence  of  the  arrival). 

For  example,  subarrays  D2  and  El  were  low  for  event  9  (Colombia)  but  not 
event  12  (Peru),  which  had  about  the  same  azimuth  but  a  larger  epicentral 
distance.  A  possible  explanation  for  this  "tuning"  effect  is  that  the  crustal 
filtering  under  a  subarray  varies  with  both  event  azimuth  and  the  angle  of 


incidence  of  the  arrival. 

For  several  subarrays  with  low  correlation  coefficients,  11-pt 
(1-sec)  Levinson  equalization  filters  were  designed  using  the  LASA  sum  as 
the  reference  trace.  A  gate  length  of  70  pts  (7  sec)  was  used.  Table  10 
lists  the  correlation  coefficients  before  and  after  equalization  and  shows 
that  a  considerable  improvement  has  been  made.  Figure  12  shows  4  sub- 
arrays  of  Event  102  (E.  Kazakh)  before  and  after  equalization.  Thus,  it 
appears  that  for  the  few  subarrays  with  anomalously  low  coefficients,  signals 
can  be  equalized  using  short  Levinson  filters. 
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Figure  9.  Event  27  (Tadzhik)  as  Recorded 
by  Subarray  Outputs 
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Figure  12.  Event  105  (G.  Kazakh),  Four  Subarray.  Beiore 
and  After  Equalization 


* 
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Amplitude  variations  of  subarray  outputs  were  generally 
quite  large  (Table  11).  Thus,  amplitude  equalization  for  all  events  was 
necessary  (and,  except  for  the  few  subarrays  with  anomalously  low  coeffi¬ 
cients,  probably  sufficient).  As  expected,  the  subarrays  having  the  largest 
and  smallest  amplitudes  varied  with  different  events. 

Table  12  lists  the  correlation  coefficient  means  and  variances, 
the  signal  degradation,  and  the  signal -to -noise  ratios  for  the  LAS  A  sum. 

Figure  13  is  a  j»lot  of  average  coefficient  versus  signal -to -noise  ratio.  No 
significant  trend  is  observed,  indicating  that  the  signals  vsed  were  large 
enough  to  prevent  the  ambient  noise  from  influencing  the  correlation 
coefficients. 

Figure  14  is  a  plot  of  average  coefficient  versus  signal 
degradation.  The  trend  observed  previously  is  not  evident  here,  which  is 
rather  surprising  because  the  two  measurements  are  roughly  equivalent. 

No  explanation  can  be  offered  for  this  discrepancy;  however,  the  signal 
degradation  across  LASA  was  small  (less  than  2  db  for  all  events). 

As  stated  in  Section  III -A,  because  of  the  excellent  waveform 
similarity  within  subarrays,  correlation -coefficient  values  for  single  seismom¬ 
eters  between  subarrays  are  similar  to  those  for  subarray  outputs.  Table  13 
compares  the  two  sets  of  coefficients  for  the  two  events  common  to  both 
studies.  The  same  gates  were  used  to  compute  both  sets.  It  can  be  seen 
that  good  agreement  exists —  both  low  and  high  values  correspond. 
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Table  11 


MAXIMUM  AMPLITUDE  VARIATION  ACROSS  ARRAY 


Maximum 
Amp  Variation 
<db) 


17.  1 
23.  1 

18.5 

17.6 

14.  1 
12.2 
11.  1 

5.6 
12.  7 

15.  0 
16.0 
17.5 
14.3 

5.  1 
19.  3 

7.6 
16.7 
14.  3 
14.  0 
21.2 


Maximum 

Subarray 


Minimum 

Subarray 
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S/N  RATIO  ON  LASA  SUM 


AVERAGE  CORRELATION  COEFFICIENT 

Figure  13.  Average  Correlation  Coefficient 
Versus  Signal -to -Noise  Ratios  On 
The  LASA  Sums  for  Events  Used 
for  Subarray  Output  Processing 
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SECTION  IV 
CONCLUSIONS 


This  study  showed  that 

•  Within  subarrays,  waveform  duplication  was  excellent. 
Variations  in  amplitude  were  sufficiently  large  to  require 
amplitude  equalization  prior  to  multichannel  processing, 
but  no  more  sophisticated  equalisation  technique  was 
necessary. 

•  Between  subarrays,  waveform  duplication  was  generally 
very  good.  Again,  amplitude  equalization  was  required 
and,  for  most  purposes,  wa<»  probably  sufficient. 

•  A  few  subarrays  had  significantly  different  waveforms  for 
some  events.  The  event  location  appeared  to  determine 
which  subarrays  exhibited  this  anomalous  behavior.  A 
possible  explanation  for  this  effect  is  that  crustal  filtering 
beneath  a  subarray  varies  with  event  location.  Short 
Levinson  filters  appeared  to  equalize  the  anomalous  wave¬ 
forms  adequately. 

•  Scattered  energy  does  not  appear  to  be  a  problem  at  the 
LASA  site. 

•  Because  of  the  excellent  intrasubarray -signal  similarity, 
the  similarity  for  single  seismometers  between  subarrays 
was  approximately  the  same  as  that  for  subarray  outputs. 
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